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Abstract. Cosmological distances as a function of redshift depend on the effective curvature 
density, Ok, via the effect on the geometrical path of photons from large scale spatial curva- 
ture and its effect on the expansion history, H(z). Cosmological time, however, depends on 
the expansion history only. Therefore, by combining distance and lookback time observations 
(or other estimates of the expansion history), it is possible to isolate the geometrical curva- 
ture contribution and measure in a model independent way, i.e., free from assumptions 
about the energy content of the universe. 

We investigate two different approaches to accomplish this task; the differential and the 
integral approach. The differential approach requires, in addition to distances, derivatives 
of distance with respect to redshift as well as knowledge of the expansion history. The 
integral approach is based on measuring the integral of the inverse of the expansion history 
via measurements of cosmic time as derived, e.g., from galaxy ages. 

In this paper, we attempt to constrain the large scale curvature of the Universe using 
distances obtained from observations of Type la supernovae together with inferred ages of 
passively evolving galaxies and Hubble parameter estimates from the large scale clustering of 
galaxies. Current data are consistent with zero spatial curvature, although the uncertainty 
of fix is of order unity. Future data sets with on the order of thousands of Type la super- 
novae distances and galaxy ages will allow us to constrain the spatial curvature fix with an 
uncertainty of < 0.1 at the 95% confidence level. 
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1 Introduction 

Current observational evidence is consistent with the universe being made up of ~ 5% bary- 
onic matter, ~ 25% non-baryonic dark matter, and ~ 70% dark energy in the form of vacuum 
energy, or equivalently, a cosmological constant [1, 2]. In this concordance model, the total 
energy density is very close to the critical energy density needed for the universe to have zero 
spatial curvature on large scales. Considerable effort is put into probing the expansion his- 
tory of the universe in greater detail in order to pinpoint the properties of the dominant dark 
energy component. Also, there are many ongoing experiments aiming at observing the dark 
matter at particle level. Less energy is being put into investigating the large scale geometry 
of the universe. 
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Most inflationary models predict that the spatial curvature should be on the order of 
|Ok| ^ 10 -5 (though there are models that allow for larger curvature) [1, 3]. As such, probing 
geometry at the largest scales can probe fundamental physics. 

Traditionally, the angular size of the fluctuations in the cosmic microwave background 
(CMB) has been regarded as the most effective way to constrain spatial curvature. The 
common explanation is that for a fixed physical size of the fluctuations, when compared to 
the case of a flat geometry, positive curvature (i.e., Ok < 0) will cause the observed angle 
to be larger and negative curvature (Ok > 0) smaller. However, since the relation between 
the physical scale and the observed angle is given by the angular diameter distance, not only 
spatial geometry but also the expansion history will affect the observed angle. Specifically, the 
reason why the observed angle is not as sensitive to the matter density (as expected for a fixed 
physical scale) is that the physical scale has a very similar dependence on the matter density 
as the angular diameter distance and thus will be factored out when calculating the observed 
angle. However, a non zero Ok affects both the spatial curvature and the expansion rate. 
Since some of the sensitivity of the CMB observations is due to the effect on the expansion 
rate, dark energy properties may also affect the observed angular scale of the fluctuations. 
Conversely, if the spatial curvature of the Universe is non zero, it can mimic evolving dark 
energy [4-7]. Because of this degeneracy, measurements of the curvature parameter, Ok, 
usually involve a parameterised description of the dark energy. Recent measurements of 
Ok, utilising combinations of different cosmological probes and simple parameterisations of 
the dark energy equation of state, indicate that the Universe is nearly flat [8-16]. However, 
some parameterisations allow large deviations from a flat Universe [17]. The relation between 
expansion and curvature is also evident from how current bounds on the acceleration depends 
on curvature [18]. 

The observed scales of CMB fluctuations are sensitive to the combination HqQk, where 
Hq is the Hubble constant. Assuming that the dark energy is in the form of a cosmological 
constant and incorporating Hubble Space Telescope (HST) priors on Hq, Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) 5-year data [1] gives -0.2851 < K < 0.0099 at 95% 
confidence level (CL). If data from baryon acoustic oscillations (BAO) and Type la super- 
novae (SNe la) are also included, the limit improves to —0.0181 < Ok < 0.0071. With the lat- 
est WMAP 7- year data [16], these limits change only marginally to —0.0133 < Ok < 0.0084. 

Barenboim et al. [8] investigate how future BAO and CMB data can be used to resolve 
the degeneracy between curvature and the dark energy equation of state, within a parame- 
terised model for w{z). In [3], it is shown how Ok can be measured with high precision using 
future CMB data combined with distance observations out to redshifts where the energy 
content of the universe is dominated by matter. The distance measures can break the degen- 
eracy between Ok and Om, and are less sensitive to the late time dark energy properties. It 
was shown in [19] how future observations of the growth of structure in the universe could 
be used to constrain the expansion history and thus, together with distance measurements 
from CMB and SNe la, can constrain Ok to sub-percent accuracy. A method to probe the 
purely geometrical part of Ok using weak gravitational lensing and BAO was devised in [20] . 
Note that, in contrast to the methods described in [3, 8, 19], this technique is completely 
independent of the expansion history of the universe and thus to the dark matter and dark 
energy densities and properties. This is also true for the methods described in this paper, 
utilising data at lower redshifts. As such, they earn their merit as independent probes of 
spatial curvature even though they will not be able to compete with the (parameterised) 
constraints on Ok from CMB observations. 
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In section 2 we outline the expansion history independent methods employed in this 
paper. These methods rely upon observable quantities introduced in section 3. We use 
the available data, described in section 4, to obtain constraints presented in section 5. In 
section 6, future constraints on Ok are investigated using simulated data sets. The results 
are discussed and summarised in section 7. 



2 Disentangling cosmic curvature and expansion history 

In a Friedmann-Lemaitre-Robertson- Walker (FLRW) universe, the expansion history at late 
times when the contribution from radiation can be neglected, is given by 

e{z) = ^ = [o M (i + zf + o K (i + zf + n DE /(*)] 1/2 , (2.1) 

where IJm, Ok and Ode are the current (effective) energy densities of matter, spatial curvature 
and dark energy, respectively, in units of the current critical density. In terms of the equation 
of state parameter, w(z), the function governing the evolution of the dark energy with time 
is 

1 + w(x) 



The parameterisation of the expansion history plays, however, no role in our method. Let us 
now define the quantity 

/w= IiS)- (2 - 3) 

which is proportional to the proper distance. We will make use of the -ffo-independent 
comoving distance 

D(z) = -1= sin \^h_I{z)] , (2.4) 
which is related to luminosity and angular diameter distances, c?l and d\, through 

D(z) = ^0^- = H d A (z)(l + z). (2.5) 

We note that distances depend on the spatial curvature, Ok in two ways; it affects the 
expansion rate of the Universe as well as the geometrical path of photons. Any effect on the 
expansion can be mimicked by a dark energy component with equation of state w = —1/3, 
corresponding to f(z) = (1 + z) 2 . We therefore seek to isolate the purely geometrical effect 
of curvature. 

The importance of not relying on any assumptions about the expansion history can be 
appreciated by noting that we do not need very elaborate dark energy models to mimic the 
effects from spatial curvature on distances out to moderate redshifts. There is an almost 
perfect degeneracy in distances already between a constant dark energy equation of state 
parameter, and the curvature, Ok- In figure 1, we show how the observed magnitude of a 
source changes with redshift and cosmological parameters. We have used a parameterisation 
of the dark energy equation of state given by 

w(z) = w + wi— — . (2.6) 

1 ~T" Z 
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Figure 1. Left panel: The sensitivity of the observed SN la magnitudes to different cosmological 
parameters 9 = {£1m, wo, wi, fix} as given by the derivatives dm/dwo(z), dm/dVl^z) etc. The 
cosmology is 17m = 0.3, f^K = 0,wo = — ljitfi = 0. Right panel: Contours at the 68.3 and 95% CL in 
the [wq, r2x]-plane from SN la observations. This is for the Union2 data set described in section 4.1 
assuming a fixed matter density Qm = 0.3 and w\ = 0. 



Since the shapes of dm/dwQ and dm/dfl-K are very similar, confidence contours from the 
redshift-distance relation for SNe la (here we use the Union2 compilation of SN data described 
in section 4.1) have a very degenerate structure when fitting for Ok and wq, even for fixed 
values of and w\. In a similar analysis using SDSS-II SN data [2], a non zero curvature 
is preferred when using the MLCS2k2 light curve fitter and assuming that dark energy is in 
the form of a cosmological constant [21]. 

Note that we assume that the spatial curvature is constant throughout space. This 
would not be the case in, e.g., so called Lemaitre-Tolman-Bondi models. However, it can be 
shown that in such models, when combining SN la data and CMB data, the observer need 
to be uncomfortably close to the center, thus violating the Copernican principle [22]. We 
also assume isotropy, i.e., that the universe looks the same in all directions since neither the 
CMB nor SNe la show any deviations from this simple picture [23]. 

From eq. (2.4), it is clear that (luminosity and angular diameter) distances are degen- 
erate with respect to curvature and expansion history. This degeneracy can be broken if the 
information from distance measurements is supplemented by the expansion history, E(z), or 
the integral, I(z). Two different approaches to measure fix can therefore be pursued. 



2.1 Differential approach 

The differential approach, originally outlined in [6], is based on differentiating eq. (2.4). We 
can then write 

_ {ED'f - 1 

Uk ~ Jyi ' y 1 - 1 ' 

where primes denote derivatives with respect to redshift. This is the curvature contribution 
to the geometry only; the curvature contribution to the expansion velocity is divided out by 
E(z). The differential approach requires not only distances, but also derivatives of distances 
as well as independent estimates of E(z). 
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Since Ok is constant for any FLRW model, eq. (2.7) measured at different redshifts 
would also offer a test of the Copernican principle [24]. Similar consistency relations for the 
concordance model, i.e. a flat cosmological constant dominated universe, can be formed by 
combining E(z), D(z) and its first and second derivatives [25]. 

2.2 Integral approach 

Measurements of the integral I(z) allows the integral approach to be followed. To gain 
insight into the relation between cosmic curvature, distance, and I(z) we can study a Taylor 
expansion of eq. (2.4) around a flat universe 



d ~ / + ^o K / 3 + y^/ 5 + o(nl). 



{21 



The distance D is to first order equal to I and Ok determines the magnitude of the higher 
order terms. Neglecting terms of order three and higher in eq. (2.8) leads to a second order 
equation with the solution 



10 



^W 1 + T5 



(2.9) 



The accuracy of this approximation decreases for increasing redshift. At z = 0.5, 1.0, and 
1.5 the error for |Ok| < 0.3 is smaller than 1, 2, and 6 parts in 10 5 , respectively. 

3 Observable quantities 

In the previous section we described how the differential and integral approach can be used 
to measure the large scale curvature of the Universe in a model independent way. We now 
discuss how D(z), D'(z), E(z), and I(z) can be obtained from observable quantities. 

3.1 Measuring D(z) and D'(z) 

The apparent magnitude, m, of a SN la with absolute magnitude M is related to the distance 
D through 



m 



M + 5 log 



Mpc 



+ 25 = M + 5 log (1 + z) + 5 log D, 



where 



• MsM + 5l0g '.n^) +25 ^ 



This can be rewritten as 



(3.1) 
(3.2) 



m — M 



In 10 



[ln(l + z) + lnD]. 



Consequently, the derivative of m with respect to redshift is 

5 



in 



In 10 



i iy_ 

JTz) + 



These two relations can be inverted to yield an explicit formula for the distance, 



D 



1 

1+1 



exp 



"In 10 



(m - M) 



(3.3) 



(3.4) 



(3.5) 
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Figure 2. Left panel: Sensitivity of curvature computed using the differential approach to different 
observables 9 = {E,m — M,m'}. A high value indicates that the observed quantity, e.g. E(z) is not 
very sensitive to the curvature and we expect the error of f2x to be large. Right panel: Sensitivity of 
curvature computed using the integral approach to the observables 9 = {I, m — A4}. The underlying 
cosmology was for both plots assumed to be J7m = 0.3 and fl\ = 0.7. 



and its derivative, 



D' = D 



In 10 



-m 



1 



5 (1 + z) 

in terms of observable quantities. We can now rewrite eq. (2.7) as 



(3.6) 



Ok = 



In 10 



1 



-m 



(1 + *) 



(1 + z) 2 exp [-0.41nl0(m- M)], 



(3.7) 



i.e., Ok = &k(E, m — A4,m'). If we can measure E, m — Ai, and m', we can thus constrain 
the spatial curvature. The error on Ok is given by 



<90k 
~~dE 



+ 



d(m — M.) 



Om-M + 



K 



dm 1 



(3i 



In the left panel in figure 2 we show how the sensitivity of Ok to the observables, = 
{E,m — _M,m/}, varies with redshift, assuming a flat cosmological constant universe with 
Om = 0.3. We expect uncertainties in E to dominate at low redshifts and uncertainties in 
m! to dominate at z > 1. Note that, since Ok is constant, we can combine data at different 
redshifts when constraining the global value of Ok- 



3.2 Measuring E(z) 

One way to measure E(z) is to use the evolution of lookback time with redshift. The -ffo- 
independent lookback time to an object at redshift z is given by 



t(z) = H (t - t z 



dx 



o (l + x)E(x)- 



(3.9) 
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By differentiating eq. (3.9) and rearranging terms, it is clear that the expansion history can 
be obtained by differentiation of the lookback time [26], 

1 dz 



E(z) 



(3.10) 



(1 + z) dr 

In appendix A, the relationship between lookback time and galaxy age is discussed. According 
to eq. (A. 4), the following equality holds 



dr 

dz 



-Hn 



dz 



(3.11) 



implying that the evolution of galaxy age, t age (z), with redshift can be used to measure the 
expansion history. 

3.3 Measuring I{z) 

Galaxy ages can also be used to measure the integral (2.3). Taking the derivative of eq. (3.9), 
rearranging terms and then integrating, we arrive at 



i'Z rz 

I(z)= / —(l + x)dx = T(z)(l + z)- / r(x)dx, 
Jo dz J 



(3.12) 



where integration by parts was used in the last step. Substituting galaxy ages for lookback 
times in eq. (3.12) comes at a price; an unknown constant T enters in eq. (A. 4), relating 
galaxy age and lookback time. 

In order to compute the integral in eq. (3.12), we need to measure i age as a function of 
redshift. If measurements below z m j n are lacking, we have to take the integral J* min ta,ge{x)dx 
into account. Fortunately this constant can be added to T. The integral I is therefore related 
to galaxy ages via 



I(z) = H 

where we have defined a new constant 

X = T + 



1- *agc(z)(l + Z) + / t age (x)dx 



^age 



(3.13) 



(3.14) 



which has to be determined from data. 

In terms of observable quantities, eq. (2.9) can be rewritten as 



n 



K 




exp [^(m-M)] 
OTzjl 



(3.15) 



i.e., J7k = £lK{h m — M). The uncertainty in a measurement of the curvature using the 
integral approach is given by 



<?n K 



( dtt K 



0-J + 



on* 



d(m - M) 



&m—M 



(3.16) 



The right panel in figure 2 shows how the sensitivity of Ok to the variables 6 = {I, m — Ai} 
varies with redshift. The error decreases with redshift and is dominated by the error in I(z). 
Figure 2 shows that for a given error in m — M, the uncertainty in Ok is larger for the 
integral than for the differential approach. 
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3.4 Differentiation of data 



The differential method relies on estimating derivatives of noisy data, SN la magnitudes to 
get m! and, in the case we do not have H(z) data, the derivative of galaxy ages to obtain 
E(z) (see section 4.2). 

Since, in a flat universe, m' gives a measure of H(z), SN la data has been used to show 
how H(z) increases with redshift. A common way to derive m! for this purpose is to use the 
method described in [27] for extracting the expansion history in uncorrelated redshift bins 
from SN la data [e.g., 18, 21, 28]. It turns out that this method is equivalent to the much 
simpler method of fitting straight lines to the m{z) data in redshift bins and use the slope of 
the fitted lines as estimates of the derivative (with corresponding errors). We have checked 
that fitting higher order polynomials to the data in a given bin does not affect the estimated 
slope appreciably 1 . 

Having nsN SNe la with magnitude uncertainty a m in a redshift bin of width Az gives 
an uncertainty in m! of 

a m r ~ °' m = (3.17) 

0.2886AzV«sn - 2 

This can be compared to the error in m — M given by 



Cfm-M = + <T 2 M , (3.19) 

where 

°^S)" 1/2 (S- ^ 

and a_M is given by an equivalent formula where nsN refers to the number of low redshift 
SNe la. 

Estimating H(z) from n age galaxy age estimates with fractional uncertainty / agc in a 
redshift bin of width Az, we obtain 

^ 0.031(1 +z) (^Y 1/2 (^.y 1 (im) . (3.21) 

H V ; V 500/ vo-5 y V o.i y 

Besides the error in H(z), we must also consider the error in the Hubble constant when 
estimating the error in E(z), 



"=*V(£) + (t>l- (3 - 22) 

When combining results from different redshift bins, one must take into account the corre- 
lations between the bins since any error in A4 and H$ will cause a systematic shift in the 
derived value of Qk- 



1 Note that this equivalent to saying that when estimating the differential using a Savitsky-Golay filter of 
degree n, the results is insensitive to n. 
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3.5 Integration of data 



Let us assume that we have measured a set of galaxy ages and their redshifts. The following 
formula can then be used to compute the integral of t age (z) numerically 



i-Zk i ^ 

F(z k ) = / t agc (x)dx ~ - 2j(4ge,j + ia g c,i-i)( z i ~ ( 3 - 23 ) 

^min J — 1 

where > z,- ii k > j and we assume zo = ^min- Values of F{zk) at different redshifts 
computed using eq. (3.23) are correlated. In appendix B we compute the covariance matrix, 
Ujk, for a model where the data points are equidistant in redshift, i.e. the number density 
is constant. We also assume that the errors are the same for all measurements. For the 
covariance matrix given by eq. (B.2), we find that Ujk < U^k, i-e., the off diagonal elements 
are comparable in magnitude to the variance a\ = Ukk- However, when computing I(z) 
using eq. (3.13) other sources of error enter 



(3.24) 



where ax is the error in X. The two last terms depends on t age and we should therefore com- 
pare the covariance matrix Ujk with the variance (1 + z) 2 cr| agc . For our model the covariance 
is inversely proportional to the number density and can therefore be beaten down by statis- 
tics. If we approximate the covariance with eq. (B.3), Ujk will be smaller than (1 + .z) 2 0"i agc 
as long as the constant number density a = dn/dz is 

Let us as an example take a = 25 corresponding to n agc ~ 40 galaxy ages evenly distributed 
in the redshift interval < z < 1.5. In that case U-u will be smaller than (1 + z) 2 a? 
by a factor 100. In the following we will therefore assume that the error in the integral 
approximation can be neglected. 



4 Data 



4.1 Type la supernovae 

We use the Union2 [29] compilation of SNe la, which is an updated version of the Union data 
set [30]. This data set contains SNe la from, e.g., the Supernova Legacy Survey, ESSENCE 
survey and HST observations. After selection cuts, the data set amounts to 557 SNe la, 
spanning a redshift range of < z < 1.4, analysed in a homogeneous fashion using the 
spectral-template-based fit method SALT2. 



4.2 Galaxy ages 

In passively evolving galaxies the stellar population was formed at high redshift and has since 
then been evolving without any further episodes of star formation. Massive galaxies located 
in high density regions of clusters have old stellar populations and are therefore believed to 
be passively evolving [31-37]. The ages of the passively evolving galaxies can be inferred 
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from their spectra using synthetic stellar population models [31, 32, 38, 39]. Spectroscopic 
dating of galaxies have been used for cosmological purposes in, e.g., [26, 40-43]. 

We use values of H(z) derived from measured galaxy ages in [43] for the differential 
approach (though look at the caveats in section 7). For the integral approach we use mea- 
surements of galaxy ages from [42] . 

4.3 Other probes of the Hubble parameter 

The Hubble parameter, H(z), can in principle be measured with other probes than galaxy 
ages, e.g., the time drift of redshifts [44] and the dipole of the luminosity distance [45]. 

In [46], a peak along the radial direction of the 2-point correlation function of LRG 
galaxies, consistent with the expected BAO signal was found. Using the BAO peak position 
as a standard ruler in the radial direction, enables a direct measurement of H{z) either at 
z = 0.24 and z = 0.43 or at z = 0.34. See, however, also [47, 48] for discussions of the validity 
of the results. 

For the Hubble constant, throughout this paper, we use Hq = 74.2 ±3.6 km s -1 Mpc -1 
as derived from the HST distance ladder observations of Cepheid variables, SNe la, and 
masers [49]. 

5 Results 

5.1 Differential approach 

Combining the Union2 SN la data set with values of H(z) derived from galaxy ages in [43] 
(again, note the caveats in section 7) and HST Hubble constant constraints, gives the result 
presented in figure 3. We have collected the data in three redshift bins: 0.1 < z < 0.4, 
0.4 < z < 0.8, and z > 0.8. Data below z = 0.1 are being used to constrain M. 

The upper panels show m — M and m! as derived from the Union2 data set. The lower 
left panel shows E{z) as derived from galaxy ages [43] with the HST value for Hq [49]. The 
dotted lines correspond to the theoretical predictions for a flat cosmological constant universe 
with f^M = 0.3. In order to put the data points at the same redshift, we have shifted each 
data point according to the expected shift given a specific cosmology where the curvature is 
a free parameter and the matter density is held fixed at SIm = 0.3. Note however that the 
results are not sensitive to this shifting, nor to the value of Om since the shift is performed 
over a relatively narrow redshift range. In each redshift bin, we perform a ^-minimisation 
to find the best-fitting value of Ok- The combined constraint on fix from all bins (taking 
the correlation between the bins properly into account, see section 3.4) is shown in the lower 
right panel, giving ft K = — 0.50±g;|? (95% CL). 

Adding H(z) derived from BAO data to the galaxy age data leads to slightly smaller 
uncertainties on the expansion history E(z) and thus also the curvature: fix = — O-SOlo'^g 
(95% CL), see figure 4. 

5.2 Integral approach 

We use galaxy ages from [42] to estimate the integral I{z) as a function of redshift. The 
relative error in t age is assumed to be 10% (L. Verde, private communication). Computing 
I(z) via eq. (3.13) requires knowledge of Hq and I. The value of I depends on the redshift 
range of the galaxy ages and when the star formation in the passively evolving galaxies 
halted, it must therefore be inferred from the data themselves. Since X has to be determined 
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Figure 3. Upper left: m — M. for the Union2 data set. Note that the error bars are smaller than 
the size of the plotting symbols. Upper right: m' for the Union2 data set. Lower left: Values of E(z) 
derived from galaxy ages in [43] and HST Hubble constant data [49]. The dotted lines correspond 
to the theoretical predictions for a fiat cosmological constant universe with J7m = 0.3. Lower right: 
X 2 - Xmin as a function of il K giving tl K = -0.50^;^ (95% CL). 



from data, we fit both Ok and I to the data using eq. (2.8) instead of computing ilx with 
eq. (3.15). The upper left panel in figure 5 shows I(z) computed for the best-fitting value of 
I assuming Hq = 74.2 km s -1 Mpc -1 . 

Since SNe la and galaxy ages are not measured at the same redshifts, the data are 
collected in redshift bins of size Az = 0.05. For the SNe la, we require at least 3 SNe per 
bin when computing the average magnitude, (m). In bins with more than one value of t age 
we use the average ((1 + z)t age — F(z)). The binning procedure results in n^n = H bins 
ranging from = 0.13 to 1.38. The upper right panel in figure 5 shows the binned values 
of m — Ai. 

Results are also sensitive to the Hubble constant. We therefore include Hq in the fit 
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Figure 4. Left: Values of E{z) derived from galaxy ages in [43], radial BAO constraints from [46], 
and HST Hubble constant data [49]. The dotted line corresponds to the theoretical prediction for 
a fiat cosmological constant universe with Om = 0.3. Right: x 2 — Xmin as a function of Ok giving 

n K = -0.5018:12 ( 95 % CL )- 



and then marginalise over it assuming a Gaussian prior to obtain constraints on Ok and I. 
Using the binned data, we fit Ok, I, and Hq to the data using the following x 2 -statistic 



where Dj = (1 + Zbinj)" 1 exp [^({rn)j - M)) and Ij =H [1- ((1 + z)t agc + F(z))j]. The 



errors . and oi i are computed using error propagation from the measurements uncertainties 
in SN la magnitudes and galaxy ages, respectively. 

The lower left panel in figure 5 shows the confidence levels in the [OK,2T]-plane a t 68.3 
and 95% CL. These constraints were obtained after marginalisation over the Hubble constant 
using the results of Riess et al. [49] as a Gaussian prior. Again we use a fixed value of M. 
determined from low redshift SNe la. 

The lower right panel in figure 5 shows X 2 ~Xmm as a function of Ok after marginalisation 
over Z. Our constraints on the curvature is Ok = 0.29!q'9 4 (95% CL). 

6 Future 

Current data clearly give quite weak constraints on Ok- In this section, we investigate the 
possibility to constrain the spatial curvature using future distance and age, or H(z), data. 
As our default future data set, we assume that we will have 2000 well observed SNe la in 
the range 0.002 < z < 1.5 with cr m = 0.1 mag. For passive galaxy ages, we assume 2000 
measurements in the range 0.1 < z < 1.5 with a precision of 0.1 Gyr. Since these numbers 
are subject to large uncertainties, we also discuss how our results scale with the total number 




(5.1) 
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Figure 5. Upper left: I{z) derived in redshift bins using galaxy ages from [42]. The best-fitting 
value of I and Hq = 74.2 km s _1 Mpc -1 were used to compute I(z) using cq. (3.13). Upper right: 
m — Ai for the same redshift bins computed using the Union2 data [29] and a value of M. fitted to 
low redshift data. Lower left: Confidence contours in the [5lK,2T]-plane at the 68.3 and 95% CL after 
marginalisation over Hq assuming a Gaussian prior. Lower right: x 2 — Xmin as a function of fly 
il K = 0.29±^J| (95% CL). 



l K giving 



and redshift distribution of the data. If not otherwise stated, we assume that Hq has an 
uncertainty of one percent. Our fiducial cosmology is a flat universe with f^M = 0.3 and 
n A = 0.7. 

6.1 Differential approach 

For our default future data set, the differential method gives the results depicted in figure 6, 
where we have Ok = 0.003 ± 0.23 (95% CL). In addition, we have assumed that we have 
on the order of 1000 SNe la at low redshifts in order to calibrate Ai, although the error 
contribution would be negligible also for a much smaller number of low redshift SNe. 

This error on is comparable to what one can get with future percent level bounds 
on the expansion rate, e.g. from radial BAO measurements, instead of the galaxy ages. In 
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Figure 6. Future constraints on with 2000 well observed SNc la in the range 0.002 < z < 1.5 
with a m — 0.1 mag and 2000 galaxy age measurements in the range 0.1 < z < 1.5 with a precision of 
0.1 Gyr. Upper left: m ~ M. as derived from SN la distances. Upper right: m! as derived from SN la 
distances. Lower left: E(z) derived from galaxy ages. The dotted lines are for a flat cosmology with 
Om = 0.3 and a cosmological constant. Lower right: The resulting \ 2 — Xmin gi vm g = 0.003 ±0.23 
(95% CL). 



[50], it is shown that a future large high-redshift spectroscopic galaxy survey covering 10 000 
square degrees could give a measure of H{z) at the few percent level in redshift bins of 
size Az = 0.2 in the redshift interval 0.5 < z < 1.5 by measuring the radial BAO signal. 
Combining these data with our default SN la data set, we typically obtain results as depicted 
in figure 7, where we have K = 0.04 ± 0.22 (95% CL). 

We note that our error scales as the inverse of the square root of the number of SNe 
and/or galaxy ages, depending on which dominates the error budget. For our default data 
set, the error contribution from E(z) and m! are comparable at low redshifts whereas at high 
redshifts, errors from m! dominate. Generally, the errors from m — M are negligible. 
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Figure 7. Future constraints on J7k with 2000 well observed SNe la in the range 0.002 < z < 1.5 
with a m =0.1 mag and constraints on H (z) at the few percent level in rcdshift bins of size Az = 0.2 
in the rcdshift interval 0.5 < z < 1.5 from future galaxy surveys [50]. Left: E(z) as derived from the 
radial BAO signal in future galaxy surveys and Hubble constant constraints. The dotted line is for 
a flat cosmology with Om = 0.3 and a cosmological constant. Lower right: The resulting \ 2 ~ Xmin 
giving Q K = 0.04 ± 0.22 (95% CL). 



6.2 Integral approach 

Applying the integral approach to our default future data set, we obtain the result depicted 
in figure 8. After marginalisation over I, we find Ok — —0.005 ± 0.054 (95% CL). Here we 
have assumed perfect knowledge of M and a 1% prior on Hq. In table 1, we show constraints 
at the 95% CL for different values of n sn and n age . From this table we learn two things: 
First, increasing only one of the number helps, but is not as effective as increasing both, i.e., 
for the numbers used, neither the errors from SNe or galaxy ages dominate completely the 
error budget. Consequently, n sn = n age is optimal. From an observational point of view it 
might however be easier and less expensive to increase one of the two numbers. Second, for 
n sn = n age , the errors scale as n -1 / 2 . This implies that reaching an error in Ok of ~ 0.01 at 
the 95% CL requires n sn = ?i agc ~ 5 • 10 . 

We have also investigated the sensitivity of our results with respect to the Hubble 
constant prior in the case where nsN = ^agc = 500. If the value of Hq was perfectly known, 
we would anticipate a measurement with an uncertainty of 0.052 at the 95% CL. Assuming a 
1% Gaussian prior increases this number to 0.098. Increasing the prior to 2% and 5% results 
in 0.110 and 0.119, respectively. Our forecast are clearly not too dependent on our choice of 
the Hq prior. 

The redshift range of the data is important. The higher the redshift, the more con- 
straining power our method have. Assuming nsN = n agc = 500 but data distributed in the 
interval 0.002 < z < 1.0 we find the uncertainty in Ok to be 0.173 at the 95% CL. Keeping 
the number of data points fixed and increasing the maximum redshift to 2 and 2.5 leads to 
an uncertainty in Ok of 0.065 and 0.054 at the 95% CL, respectively. 
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Figure 8. Future constraints on fix with 2000 well observed SNe la in the range 0.002 < z < 1.5 
with <j m ~ 0.1 mag and 2000 galaxy age measurements in the range 0.1 < z < 1.5 with a precision 
of 0.1 Gyr. Upper left: I(z) as derived from galaxy ages. Upper right: m — M. as derived from 
SN la distances. Lower left: Confidence contours in the [r2x,I]-plane at the 68.3% and 95% CL after 
marginalisation over Hq assuming a 1% Gaussian prior. Lower right: The resulting x 2 ~~ Xmin §i vm S 
n K = -0.005 ± 0.054 (95% CL. 



7 Summary and Discussion 

We have investigated how the large scale spatial curvature of the Universe can be measured 
model independently. Cosmological distances depend on Ok via the expansion history and 
the geometrical path of photons. If the two contributions can be disentangled, fix can be 
measured independently of the energy content of the Universe. We have pursued two different 
paths to disentangle the effects of geometry and expansion. 

The differential approach, originally suggested by Clarkson, Cortes, and Basset [6], 
allows Ok to be measured via eq. (2.7). In addition to distances, this approach requires 
derivatives of distance with respect to redshift and independent measurements of the ex- 
pansion history. The expansion history can be measured using different probes. Here we 
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Table 1. Constraints on fix at the 95% CL from simulated SN la and galaxy age data. The data is 
uniformly distributed in the redshift interval 0.002 < z < 1.5. For SNe la and galaxy ages the error 
is a m = 0.1 mag and er ta e = 0.1 Gyr, respectively. The value of the Hubble constant was assumed to 
be known to 1%. 





n age = 100 


n age = 500 


n agc = 1000 


n agc = 2000 


100 


0.188 


0.188 


0.165 


0.156 


500 


0.122 


0.098 


0.093 


0.084 


1000 


0.111 


0.080 


0.072 


0.066 


2000 


0.104 


0.069 


0.060 


0.054 



focus on the method proposed by Jimenez and Loeb [26], which utilises relative galaxy ages. 
Cosmological distances with corresponding redshift derivatives can be derived from SN la 
observations. 

The differential approach requires numerical derivation of data to obtain D'(z), and in 
the case of galaxy ages, also E(z). Numerical derivation of noisy data is far from trivial and 
we suspect that errors have often been underestimated in the past. In [42], for example, they 
use a customised method to obtain H(z) that gives uncomfortably small errors. Basically, 
they first group together all galaxies that are within Az = 0.03 of each other to get an 
estimate of the age of the universe at a given redshift with as many galaxies as possible. 
The redshift interval is chosen to be small enough to avoid incorporating galaxies that have 
already evolved in age, but large enough to have more than one galaxy in most of the bins. 
Age differences are computed only for those bins in redshift that are separated by more than 
z = 0.1 but not by more than z = 0.15. The first limit is imposed so that the age evolution 
between the two bins is larger than the error in the age determination. The authors claim 
that this procedure provides a robust determination of dz/dt age . We have not been able 
to reproduce their results and have therefore not been able to check the stability of this 
claim. It is, however, unlikely that the given uncertainties in the derived values of H(z) are 
correct. From a sample of 32 galaxy ages (out of which some are discarded as outliers), you 
could form a maximum of 16 pairs of 2 galaxy ages each which could give 8 values of the 
Hubble parameter. Assuming an error of 10% in the galaxy ages, a simple analytical error 
propagation gives an error of > 100% in H(z). A straight line fit in bins give similar results; 
at best a 100% error on H(z) at five different redshifts. With this (lack of) precision it turns 
out to be futile to determine Qk- In [43], dt &ge /dz is computed by defining the edge in the 
t agc vs z-plot formed by the oldest galaxies in a sample of 24 clusters, yielding two additional 
H(z) estimates at z = 0.48 and 0.9. We have not checked the stability of this method but 
the derived H(z) uncertainties seems to be more reasonable when compared to the number 
of galaxy ages employed. 

Galaxy ages as a function of redshift can, alternatively, be used to measure the integral 
I(z) via eq. (3.13). Knowledge of I(z) allows f^K to be measured via the integral approach, 
see eq. (2.9). This approach requires numerical integration of noisy data instead of numerical 
differentiation. There are potential problems also with the integral method. The integral 
approach, which relies on measurements of lookback time, is not as general as the differential 
approach, where, in principle, E(z) can be measured using several fundamentally different 
methods. We assume that the lookback time can be measured via passively evolving galaxies 
acting as cosmic chronometers. When converting galaxy ages to lookback times we make the 
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assumption that the star formation halted at the same time, i.e., that all passively evolving 
galaxies share the same birthday. Even if a population of cosmic chronometers exists, there 
is always a risk that our sample is contaminated by other galaxies. Two unknown quantities 
enters our calculation of I(z): the constant offset in time between galaxy age and lookback 
time and the low redshift part of the integral over galaxy age. Fortunately, these unknowns 
can be handled by a single nuisance parameter I. The need for an extra nuisance parameter 
is however another drawback of the method as it degrades the constraining power of the data. 
Another problem with the integral method is that values of I{z) computed using eq. (3.23) 
are correlated. However, as shown in appendix B the correlations decrease with increasing 
number of galaxy ages. 

Generally, the differential method is the most effective if we have, in addition to distance 
measures, independent H(z) constraints from, e.g., radial BAO measurements. The integral 
method, on the other hand, outperforms the differential method if the expansion history 
needs to be reconstructed from lookback time measurements. The reason why this is not 
the case for current data is that the derived H(z) from passive galaxy ages probably have 
underestimated uncertainties. 

The model independent constraints on fix we present here are rather weak, — 1 < fix ^ 
1 at the 95% CL, and at least for the differential method, there are reasonable doubts that 
the uncertainty in the derived expansion history has been underestimated. We therefore 
investigated to what degree a proper analysis of future data could sharpen current curvature 
constraints. Even if the amount of data is significantly increased, we find that the uncertainty 
in fix w iH be fairly large, ~ 0.1 at the 95% CL, at least compared to (model dependent) 
CMB constraints. It is therefore unrealistic that the methods presented in this paper will 
be able to constrain, e.g. inflationary models that predict spatial curvature on the order of 

I^kI < io- 5 . 

In the introduction we noted another proposed method of measuring fix model inde- 
pendently using weak gravitational lensing and BAO by Bernstein [20]. Forecasts for this 
method, assuming a full sky survey of galaxies in the redshift range < z < 3 supplemented 
with photometric redshifts, predicts an uncertainty in fix of ~ 0.04 at the 68.3% CL. This 
precision is comparable to what could be obtained with the integral method with n sn = 1000 
and n age = 500 (see table 1). 
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A Lookback time and galaxy ages 

The integral I(z) can, according to eq. (3.12), be obtained from measurements of lookback 
time as function of redshift. We assume that the lookback time can be measured via the age 
of galaxies, or rather the age of the stellar population in galaxies. 
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To understand the relationship between lookback time and galaxy age, we write the 
present age of the Universe as 



t -t oo = - —dz (A.l) 
= (t - t g ) + (t x - i c ) + (t c - too), (A.2) 

where to — t z is the lookback time, t z — t c is the age of a galaxy at redshift z, which was 
created at redshift z c , and t c — is the age of the Universe when the passively evolving 
stellar population of the galaxy was created. 

We naively assume that all the galaxies share a common time of star formation. In that 

case 

(t -t z ) = T-{t z -t c ), (A.3) 

where the constant T is the age of the Universe at redshift z c subtracted from the present 
age of the Universe. The lookback time and the galaxy age are hence related by an additive 
constant. Consequently galaxy ages are related to r through 

T = tf (T-tage). (A.4) 



B Derivation of error in the integral approximation 

We use eq. (3.23) to compute the area under the graph of t a _ ge (z). Let us make some sim- 
plifying assumptions that will allow us to draw some conclusions about the accuracy of the 
approximation. We assume all the measurements i a ge,i at redshifts z% to be separated by 
Az. Furthermore we assume the error o"j ago to be the same for all £ a ge,i- Values of F(z) 
computed at different redshifts using eq. (3.23) are obviously correlated. With the previous 
assumptions, we can derive an approximate formula for the covariance matrix using error 
propagation, 

AzV? < min(?', k)Az 2 a? . 

(B.l) 

The assumption that the points are evenly distributed in redshift means that dn/dz = a is 
a constant. Consequently k = a(zk — z m \ a ) and Az = a" 1 . We can therefore rewrite the 
covariance in terms of redshift as 

Ujk < [mm(zj,Zk) - z min ] . (B.2) 
At redshift Zk the variance is hence 

4 = U kk <^(z k - z min ) . (B.3) 
a 

We note that Ujk ^ U kk - 
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